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A recent discovery concerning the importance of hadron production in space radiation 
is that pions can contribute up to twenty percent of the dose from galactic cosmic ray 
interactions (S. Aghara, S. Blattnig, J. Norbury, R. Singleterry, Nuclear Instruments and 
Methods, Vol. 267, 2009, p. 1115). Although the contribution for dose equivalent will 
be smaller, the dose contribution could be important for fluence based radiation models. 
Pion production cross sections will be an essential ingredient to such models, and it is 
of interest to investigate the adequacy of the pion production experimental data base for 
energies relevant to space radiation. The pion production threshold in nucleon - nucleon 
reactions is at 280 MeV and, in an interesting accident of nature, this lies near the peak of 
the galactic cosmic ray proton spectrum. Therefore, pion production data are needed from 
threshold up to energies around 50 GeV/nucleon, where the galactic cosmic ray fluence is 
of decreasing importance. Total and differential cross section data for pion production in 
this energy range will be reviewed. The availability and accuracy of theoretical models will 
also be discussed. It will be shown that there are a significant lack of data in this important 
energy range and that theoretical models still need improvement. 


Nomenclature 


p Proton 

n Neutron 

d Deuteron 

7r Pion 

p Muon 

v Neutrino 

e Electron 

7 Photon 

r Lifetime, sec 

c Speed of light, m/sec 

A Nucleus, or Nucleon number 

M Meson 

X Any particle 

s Square of center of momentum energy, GeV 1 2 

T Kinetic energy, GeV 

E Total energy, GeV 

p Momentum, GeV (when c = 1) 

Pt Transverse momentum, GeV (when c = 1) 
a Total cross section, mb 

sy Lorentz invariant differential cross section, GG ™^ sr (when c = 1) 
p Density, kg/m 3 

N a Avagadro’s number 

l Mean free path 

* Senior Research Scientist and Leader of the Space Radiation Group; Durability, Damage Tolerance and Reliability Branch, 
NASA Langley Research Center, Hampton, Virginia, 23681, USA, AIAA Senior Member. 
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I. Introduction 


The space radiation environment comprises energetic particles produced from three sources. Each source 
has a characteristic spectrum, which is the number of particles as a function of energy. Firstly, solar particle 
events (SPE) consist primarily of protons emitted from the Sun during coronal mass ejections and solar 
flares. These events are rare, but when they occur they can potentially inflict a lethal dose of radiation to 
astronauts. They are also of great concern for the stability of electronic devices. Energies can reach up to the 
GeV region. Secondly, galactic cosmic rays are protons and heavier nuclei emitted from supernovae within 
our Milky Way Galaxy and accelerated to the vicinity of the solar system and elsewhere. The particles have 
energies from 10 AMeV out to the ZeV region. (Zetta eV = 10 21 eV.) The peak of the spectrum is in the 
MeV - GeV region, which is most important for space radiation. For space radiation purposes, it is sufficient 
to consider the spectrum out to only about 50 GeV, because the flux is so much smaller beyond that energy. 
Thirdly, the geomagnetically trapped particles are protons and electrons confined by the magnetic field of 
Earth. There are two distinct regions called the inner and outer Van Allen radiation belts. Protons and 
electrons are found in both belts, with the most important being inner belt protons with energies up to 
100 MeV and beyond, and outer belt electrons with energies up to 100 keV and beyond. It is important to 
have an accurate knowledge of the spectrum of these particles at various places in the solar system, where a 
particular mission will take place. The spectrum will be modified when the particles impinge on a medium. 
For instance, the spectrum of particles will be modified after traversal through a spacecraft wall, due to 
the many atomic and nuclear reactions that take place. Knowledge of the spectrum on the other side of 
the medium enables one to determine the dose of radiation that an astronaut receives. Several computer 
codes have been developed at NASA Langley Research Center that enable one to calculate the spectrum 
of particles on the other side of a medium, and the subsequent radiation dose. Codes at other institutions 
have also been developed, but they are often limited to a small number of particles and also, due to their 
Monte Carlo nature, take a long time to run. An important feature of the Langley codes is that they are 
also able to deal with the spectrum of heavy nuclei found in the galactic cosmic rays. A unique feature of 
the Langley codes is their short run time and ease of use, making them ideal for engineering design. The 
codes are deterministic and are able to run quickly allowing for almost real time simulations. Nuclei with 
low nucleon number provide the best protection against space radiations. This means that, for example, 
liquid hydrogen, water and polyethylene provide much better protection than materials such as aluminum or 
lead. This has major implications for human space missions and provides insight as to how to best protect 
astronauts and electronic devices. Transport codes are also widely used in studies of cosmic rays propagating 
through various media. If one knows the cosmic ray spectrum incident on top of the Earth’s atmosphere, then 
one can deduce the cosmic ray spectrum observed on the ground by transporting through the atmosphere. 
Similarly, if one knows the cosmic ray spectrum incident on a spacecraft wall then one can deduce the 
radiation environment inside a spacecraft. 1 Deterministic transport codes are concerned with solving the 
Boltzmann transport equation, whereas Monte Carlo codes simulate particle trajectories. HZETRN (High 
Z and Energy TRaNsport) is a computer code that is able to calculate the radiation environment inside a 
spacecraft, given knowledge of the external environment and the reaction cross sections. It is a deterministic 
space radiation transport code based upon the straight-ahead approximation. 1 

II. Hadron production &; interaction 

Hadrons are produced in proton-proton, proton-nucleus and nucleus-nucleus collisions. In a high density 
shield, such as a spacecraft wall, the hadrons are typically absorbed before they decay. In a low density 
environment, such as the Earth or Mars atmosphere, the produced hadrons decay before further interactions, 
leading to a buildup of muons, electrons and photons, which is called an electromagnetic cascade. It has 
recently been found that pions can contribute up to twenty percent of the dose from galactic cosmic rays. 2 
(The dose equivalent has not yet been considered.) The kinetic energy threshold for pion production in 
proton-proton collisions is 280 MeV, and from the energies discussed previously, one can conclude that 
hadron production will be important for cosmic rays and possibly some energetic solar particle events, but 
not for geomagnetically trapped particles. Note that pp cross sections are quite large. The cross section is 
about 50 mb in the intermediate energy region and grows much larger as energy is increased. 

It is now over seventy years since the prediction of the pion by Yukawa, and since then hundreds of 
other hadrons have been discovered and studied extensively. The pion is the most important hadron because 
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it is the lightest, and therefore most easily produced and is responsible for the longest range part of the 
nucleon-nucleon interaction. One might imagine that there are no major gaps remaining in the experimental 
study of the pion. However, this is far from the case. The region of maximum galactic cosmic ray intensity 
is the 0.1 - 10 GeV region, yet the energy region with the least amount of data are in the 1-10 GeV range. 
This is called the intermediate energy region. Furthermore, the region where theory is least understood is 
1-10 GeV because intermediate energy is too high for non-relativistic models to work well and too low 
for ultrarelativistic approximations. This is made clear in the work of Eyser and Machleidt, 3 who analyze 
spin correlation functions and show that the standard one boson exchange models, fail badly in explaining 
spin correlation cross section measurements. Nevertheless, they show that these models do work reasonably 
well when spin averages are taken into account. For space radiations studies, of course one does not need 
spin cross sections, but the point is that 70 years after Yukawa, there are still problems with interpreting 
intermediate energy measurements. 

Examples of exclusive cross sections are 

p + p — > 7r + + p + n, (1) 

p + p — > 7r + + d, (2) 

p + p — ► 7T + + p V n V 7T°, (3) 

where all final state particles are identified. In an inclusive cross section, these are all summed over as in 

p + p — > 7r + + X, (4) 

where X represents anything. Calculating an exclusive cross section is easier than an inclusive one, because 
one has to add several or many exclusive calculations to obtain the final exclusive cross section. Experimen- 
tally, it is much easier to measure only one particle, such as the 7r + , rather than many particles. 

The FLUKA and MCNPX codes both include meson production. The FLUKA code uses different 
interaction models depending on the energy region; below 2.5 GeV the PEANUT (PreEquilibrium Approach 
to Nuclear Thermalization) model including a Generalized IntraNuclear Cascade, from 2.5 to 4 GeV a 
resonance production and decay model including a Glauber Cascade, above 4 GeV the Dual Parton Model. 
The DPM is based on the more general Gribov-Regge theory (GRT) formulated in terms of string interactions. 
The Gribov-Regge theory can be implemented in a variety of ways, which give rise to the different models 
such as the DPM, the Very Energetic Nuclear Scattering (VENUS) model or the Quark Gluon Strong Model 
(QGSM). The Los Alamos version of the quark gluon string model is called LAQGSM. The AMOR (Artru- 
Mennessier Off-Shell Resonance) model is used in VENUS, but other possibilities are the Field-Feynman 
model or the Lund JETSET model. Others include multi-peripheral, parton cascade and string-parton 
models. One problem with many nuclear interaction models is that they use Monte Carlo procedures leading 
to computational inefficiency. Parametric methods can relieve the computational burden. Meson production 
is now included in a research version of HZETRN with pion and muon transport. The plan is to include 
neutral pion and eta meson production, which are the main source of the electromagnetic showers and to 
include charged kaons which are additional sources of muons. 

Particle production in the high energy region is also well described by several types of arithmetic param- 
eterizations, which incorporate the concept of Feynman scaling. Their advantage is that they are extremely 
fast to compute and are therefore suitable for transport codes, such as the deterministic code HZETRN, 
which run fast enough to be useful for engineering design. Therefore, studies 4,5 have been made concerning 
the accuracy of these parameterizations for production of pions, kaons, protons and antiprotons in proton- 
proton collisions. Figure 1 shows that these can be very accurate. Also, note that these are not mere 
parameterizations of data. They are based on high energy physics models incorporating important physical 
ideas such as scaling, which comes about due to the parton nature of hadrons. These parameterizations 
are also useful in applications to proton-nucleus 4 and nucleus-nucleus high energy reactions. Figure 1 also 
shows pion production from high energy proton-Carbon collisions. Again the agreement is very good for 
some parameterizations. However, at lower energies, these parameterizations fail, especially in describing 
differential cross section measurements 5 for nucleus-nucleus collisions. This is shown in Fig. 2, which also 
shows that a simple thermal model parameterization gives excellent results. 5 

In the lower energy region, FLUKA and other codes use a boson exchange model. This model does 
not include spin and does not include all the higher meson exchanges explicitly, but rather includes pa- 
rameterizations of a simple scalar model. A model which does not include spin should really be based on 
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Figure 1. 7r+ production in pp (left panel) and pC (proton-Carbon) collisions (right panel). The values of 

the Feynman scaling variable xf from top to bottom are 0.0, 0.01, 0.02, 0.025, 0.03, 0.04, 0.05, 0.06, 0.07, 
0.075, 0.08, 0.1, 0.12, 0.15, 0.2, 0.25, 0.3, 0.35, 0.45, 0.55. Following reference 6 data and lines are multiplied 
successively by 0.5 to allow for a better separation. The parameterization has been multiplied by a best 
fit factor of 12 0,9 . Data from references 6 ’ 7 are plotted against the parameterization of Mokhov. 8 The pC 
parameterization has been scaled by the Cronin enhancement factor as discussed in reference. 4 Reproduced 
from reference 4 with permission. 


a fully consistent scalar field model and at intermediate energies one might expect such a theory to work 
well. Such a scalar field model is able to explain available data quite well for both 7 r and 77 production . 9 
The dominant reaction mechanism for 77 production is via formation and decay of the iV*(1535) resonance 
through pp — > pN* (1535) — > pprj. Future work will involve implementing a scalar field model 9 into HZETRN, 
with the higher energy parameterizations being used only above intermediate energy, in a manner similar to 
FLUKA and GEANT, which use different models at different energies. 

In the very low energy region near the pion and kaon threshold, there are interesting things to consider . 10 
Particle reaction transport codes, such as GEANT and FLUKA are widely used in the design of accelerator 
experiments and the simulation of particle detectors. While typical cross sections away from threshold reach 
tens of millibarn, the cross sections near threshold are about a thousand times smaller, typically in the 
microbarn region. The fact that theory and experiment disagreed so much in the threshold region was 
therefore of no real concern for particle simulation codes because the cross sections in the threshold region 
are so small. An interesting accident of nature is that the peak of the cosmic ray spectrum, for both protons 
and heavier nuclei, occurs near the pion production threshold. The Boltzmann transport equation contains a 
term which is the cosmic ray flux multiplied by the cross section. Therefore, when considering pion and kaon 
production from proton-proton reactions, small cross sections at low energy can be as important as larger 
cross sections at higher energy. This is also true for subthreshold kaon production in nuclear collisions, but 
not for subthreshold pion production. 
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Figure 2. Inclusive 7r — cross section for Ar+KCl collisions. Badhwar parametrization 11 (solid line) multiplied 
by ApAx fails to describe data at 800 MeV/nucleon. Lab angles are indicated. Even if the multiplication 
factor is changed to an arbitrary constant, the theoretical shape does not match the shape of the data. Other 
parametrizations are of similar poor quality. A thermal spectrum (dashed line) successfully describes data. 
Reproduced from reference 5 with permission. 
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Both the nuclear and particle reactions can be written together compactly as 


A + A 

-> A' + AT, 

(5) 


-A M + X, 

(6) 

M + A - 

-> A' + X, 

(7) 

- 

-a M + X, 

(8) 


where A is any nucleus, including a proton or neutron, and M is a meson such as 7 r + , n~ , 7 r°, K°, K + , 
or K~ . The neutral kaon contains an equal contribution from the short and long components, as in K° = 
50%A/g + 50%Ar/. Note that A' can also be A. The nuclear reactions (5) do not involve a meson M on either 
the left or right side of the reactions, whereas the particle reactions (6), (7), (8) do involve M on either the 
left or right side of the reactions. Such a compact way of writing the reactions is useful when keeping track 
of particles in a transport code. Note the vast number of possible reactions, because A can be any nucleus, 
and the pions and kaons can be charged or neutral. 

When pions or kaons are produced, they decay according to the following major decay modes. The 
charged pions undergo a weak decay 12 

7T* — > (99.99%) r = 2.6 x 1CT 8 sec, cr = 7.8 m, (9) 

where 99.99% refers to the percentage of this mode of decay relative to all the possible decay modes. The 
lifetime r has been listed as well as the quantity cr. Both of these quantities are in the rest frame of the 
decaying particle. The quantity cr gives the approximate distance the pion will travel before decaying. 
This is in the rest frame of the pion. Of course in the lab frame, which is the target frame or spacecraft 
wall frame, the lifetime will appear longer and the distance will therefore be longer. Thus, cr is actually 
a minimum distance. Given that cr ss 8 m for the above reaction, then the charged pion will interact and 
produce ionization before decaying. For this reason, one must also include the stopping powers of charged 
pions in radiation transport codes. The neutral pion undergoes an electromagnetic decay 12 

7T° — ■> 2y (99.8%) r = 8.4 x 10 _17 sec, cr = 25 nm. (10) 

Thus, the neutral pion can be considered not to propagate at all, but to instantaneously decay into two 
photons 7 . These photons will immediately decay via 

7 — > e + + e~ . ( 11 ) 

The electrons will produce bremsstrahlung and ionization. The positrons will annihilate, 

e + + e - — > 27 , ( 12 ) 

to produce more photons which will decay, producing more electrons and positrons. Thus, a cascade of e + e - 
pairs and photons will be produced from the decaying neutral pion, which can be schematically written as 

7 T° — > electromagnetic cascade = e + , e _ , 7 - (13) 

The above considerations show why it is important to include pions and kaons in radiation transport 
codes. The particle flux as a function of Earth atmospheric depth steadily decreases. For typical jumbo jet 
altitudes (30,000 feet or 10 km), the atmospheric depth is about 200 g/cm 2 , whereas at altitudes appropriate 
for supersonic flight (20 km), the atmospheric depth is only about 50 g/cm 2 . Even though neutrons are not 
present in the external cosmic ray spectrum, nevertheless, neutrons are dominant even at zero atmospheric 
depth. This is due to the neutron albedo from the atmosphere, and neutrons dominate the particle spectrum 
almost all the way down to the Earth surface. This, and the large relative biological effectiveness value, 
is why neutrons are so important for radiation received in aircraft. The reason for the monotonic fall in 
particle flux as a function of atmospheric depth is due to absorption from nuclear reactions. However, note 
that for small depths the pion and muon flux increase. Also, note that the muon flux is even larger than 
the proton flux at larger depths. Note that the large electron flux is due to 7 r° decay and the associated 
electromagnetic cascade. (Electrons due to atomic ionization have much smaller energies.) This feature of 
the buildup of electromagnetic particles (electrons, photon, muons) is called an electromagnetic cascade and 
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is important for radiation transport through less dense media such as Earth and Martian atmospheres. Now 
consider why this buildup of pions and electromagnetic particles is so important for atmospheres. It also 
may be important for dense media such as spacecraft shielding. 

Recall that one mole is the number of atoms in a 12 g sample of 12 C. This number is determined from 
experiment to be Na = 6.02 x 10 23 , which is called Avagadro’s number. The interest is in nuclear reactions 
and so the number of nuclei is set to be approximately the number of molecules. Also, assume air to be 
composed entirely of Nitrogen. A 27 g sample of 2 'A1 or a 14 g sample of 14 N contains Na molecules. Let 
p be defined as the number of nuclei per cubic centimeter, 


P = 


# nuclei 


cm J 


pN A 
A ’ 


(14) 


where p is the usual mass density (g/cm 3 ) and A 

is the atomic number. The mean free path / , which is the 

average distance between nuclear collisions, is given by 13 


l = — 

pa 

A 

(15) 

paN A ' 

The mean free path l will have units of cm if p 

is expressed in g/cm 3 and a in cm 2 . 

Typical spacecraft 

shields are made from Al, tissue can be approximated as water, the Earth atmosphere can be approximated 
as Nitrogen, and the Mars atmosphere is mostly CO 2 with a total aerial density of about 10 g/cm 2 . For 

these materials, 



PAl = 

3 g/cm 3 , 

(16) 

Pwater — 

lg/cm 3 , 

(17) 

Pair — 

10 -3 g/cm 3 , 

(18) 

PMars atmosphere 

10 _ 5 g/cm 3 = (l/100)p air . 

(19) 

There are rough approximations made above, but only order of magnitude estimates are 

sought. Define 

h = 

A 

(20) 

pa ± N A ’ 

ai = 

lmb = 10" 27 cm 2 , 

(21) 

and the following results are obtained 



l 

'l Ai « 150 m, 

(22) 

W water ~ 300 m, 

(23) 

i 

1 air « 200 km, 

(24) 

W Mars atmosphere ^ 20, 000 km. 

(25) 

At a momentum of 1 GeV, 12 



°pp 

, w 50mb, 

(26) 

^7T+p 

, « 20mb, 

(27) 

so that 



Iai(pp) 

« 150 m/50 = 3 m, 

(28) 

Lvater ( PP ) 

« 300m/50 = 6m, 

(29) 

kir (PP) 

« 200 km/50 = 4km, 

(30) 

Gla.rs atmosphere (PP) 

« 20, 000 km/50 = 400 km, 

(31) 

and 



/ A l(7T + p) 

« 150 m/20 = 7.5 m, 

(32) 

^water(?Ap) 

« 300 m/20 = 15 m, 

(33) 

Gir (7r + p) 

w 200 km/20 = 10 km, 

(34) 

Glars atmosphere (A P ) 

w 20, 000 km/20 = 1000 km. 

(35) 
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The numbers above are all rough order of magnitude estimates. Recall the decay length of cr = 7.8 m. 
Multiplying by the relativistic 7 ss 4 factor for a 1 GeV pion gives ycr ss 30 m. Thus, the conclusion is that 
in a medium such as an atmosphere, a charged pion will decay into a muon before any proton interactions 
occur, leading to a buildup of muons. However, in condensed matter (such as A1 or water), this will happen 
less often. This is the reason for the large buildup of charged pions and muons in atmospheres. Similar 
arguments can be made for neutral pions ( 7 r°), which have cr = 25 nm. For all types of media, these will 
all decay to two photons before being absorbed. Thus, for all media, 7 r° production and decay must be 
considered. For a 1 GeV pion, the photons will have energies in the GeV region. Consider the range of 
photons as a function of aerial density. GeV photons reach depths as high as 100 g/ cm 2 . However, photons 
have a relatively small biological effect compared to neutrons or heavy ions. In conclusion, neutral pion 
production and decay must be included for transport through shields and atmospheres. The issue of the 
importance of charged pions depends on the transport medium . 2 Both neutral and charged pion production 
and decay must be included for transport through atmospheres. 

III. Total cross section for single pion production 

Apart from the neutron and proton, the pion is the most important hadron because it is the lightest 
meson and is therefore the most copiously produced hadron. The most fundamental observable describing 
its production is the total cross section for producing a single pion. Given the hundreds of other hadrons 
and the dozens of other experimental observables, one would assume that this fundamental observable, the 
total cross section for producing a single pion, is thoroughly understood. It will be shown that this is far 
from the case. 

One of the most important conservation laws in strong interaction physics is isospin conservation, which 
implies the relations between cross sections as summarized in table 1. This shows that if one has experimental 
data for certain cross sections, then data for all particle cross sections need not be measured, because isospin 
conservation can be used to obtain equivalent cross sections for the remaining reactions which have not been 
measured. Similarly, if one does have data for these remaining cross sections, then the data can be used to 
test the isospin predictions. Consider the relation 

a(pn — > ppix~) = a(pp —> ppK°) — -cr(pp —> pmr + ) + tj(pn —> pmr°), (36) 

which is obtained from table 1 because a(pn — > ppix~) = a(pn — » rm7r + ). 

Table 1. Isospin relations between single pion production cross sections. If one assumes that data are provided 
for the reactions indicated, then isospin conservation provides equivalent cross sections for the remaining 
reactions. 

Reaction Equivalent cross section 

pp — > ppn° data 

pp — y pmx + data 

pn — » pmx° data 

pn — > nmx + cr{pp — > ppn 0 ) — \<j{pp — > pnix + ) + a(pn — > pmx°) 
pn — > ppn~ a(pn — > nnix + ) 

nn — > nn7r° cr(pp — > ppn 0 ) 

nn — > npn~ cr(pp — > pnix + ) 

pp — y g?7t + data 

pn — > dn 0 (pp — > dn + ) 

nn — y cbr~ u(pp — > (i7r + ) 
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In order to aid the reader, table 2 gives the threshold energy for the reaction pp — > ppir °, and table 3 
lists the relation between kinetic energy T and the root of the center of momentum totals energy y/s. 

Table 2. Threshold energy for pp — > ppn° . 

Threshold parameter Value (GeV) 


T\ a h 

0.28 

-Slab 

1.2 

P\&b 

0.78 

A 

2.0 


Table 3. Relation between T and ^ s . 

T (GeV) y/s (GeV) 

1 2.3 

2 2.7 

3 3.0 

4 3.3 


Data for cr(pp — » ppn 0 ), u (pp — > pnir + ), a(pn — » pmr° ), a(pn — > ppn~) are given by Teis et al. 14 and 
Bystricky et al. 16 Cross section values, at a variety of energies, are taken from these references and are listed 
in table 4, which clearly shows that the experimental data are not consistent with Eq. (36). This is a very 
poor state of affairs for such a fundamental observable. These isospin relations are also true for differential 
cross sections. Transport codes, such as FLUKA and GEANT, rely on such isospin relations to calculate 
particle production. Data are needed to resolve this issue. 

Table 4. Test of the isospin relation a(pn — > ppir~) = rr(pp — > pp 7r°) — \&(pp — > pn7r+) + a (pm — > pmr°), which implies 
that the numbers in the last two columns should be equal. 

y/s cr(pp — > ppir 0 ) a(pp — > pnir + ) a(pn — > pmr°) <r(pn — > ppn~) &(pp — > ppir 0 ) 

— ^tj(j>p — » pmr + ) 

+a(pn — y pmr°) 


2.1 

0.3 

2 

1 

0.3 

0.3 

2.15 

1 

5 

3 

0.5 

1.5 

2.2 

3 

12 

7 

1.5 

4 

2.3 

3.5 

17 

14 

2.2 

9 

2.4 

3.7 

18 

7 

2.5 

1.7 

2.5 

3.5 

17 

7 

2.6 

1 


IV. Summary of pion data 

Consider the status of experimental data for pion production. Even though there are hundreds of other 
particles to consider, the poor situation for pions illustrates the problem for all particles, especially as pions 
are better measured and understood compared to other particles. Table 5 summarizes the available data for 
pion production in nucleon - nucleon collisions. This will now be considered. In the discussion below, all 
energies refer to kinetic energy (i.e. T, not y/s) unless otherwise stated. 


9 of 11 


American Institute of Aeronautics and Astronautics 



A. Exclusive single pion production 

1. Total cross section 

A lot of the data are summarized in the references. 14-17 While there are quite a few measurements available, 
the following problems are noted. Many data do not go beyond a few GeV in kinetic energy. All exclusive 
channels have not been measured. Of the exclusive channels that are measured, some of the measurements 
appear inconsistent with isospin conservation, which is used in codes such as FLUKA to calculate some 
cross sections. Single pion production data, 14 which have been shown above to be inconsistent with isospin 
conservation. Also, note that much of the data are for low values of kinetic energy. There are also difficulties 
in predicting data near threshold. 18 

2. Differential cross section 

The isospin relations in table 1 also hold for differential cross sections, which need to be checked to ensure 
that they do not violate isospin constraints. There are very few data and most are below 1 GeV. Again, 
there are difficulties near threshold. 18 Although, not of interest for space radiation, Eyser and Machleidt 3 
have shown that meson exchange models are unable to predict spin correlation functions in the GeV region. 

B. Exclusive multiple pion production 

A lot of data are summarized in the references. 16, 17, 19 The experimental situation for total cross section is 
poor with few measurements available. There do not appear to be any data available for differential cross 
sections. 

C. Inclusive pion production 

Inclusive measurements are much easier than exclusive measurements, and one would think that the data 
situation is much better. However, most of the data are at very high energy. In the intermediate energy 
region (1 - 10 GeV), there are very few measurements of either total or differential cross sections. This is 
not a good situation for space radiation, since it is the inclusive cross sections in the intermediate energy 
region that are most often required. One is therefore forced to use high energy parameterizations and hope 
that they work at intermediate energy. Note that these problems in the intermediate energy region are not 
of much concern for many transport codes, such as GEANT and FLUKA, when used for accelerator and 
detector design for particle physics experiments, because the intermediate energy region is small, and there 
are many other energy regions which are just as important. However, for space radiation purposes, where the 
cosmic ray spectrum peaks at intermediate energy, then these problems are very important. The discussion 
has focused on pion production and the situation is worse for the other hadrons. 


Table 5. Summary of pion production data from nucleon - nucleon collisions. Single, multiple and inclusive 
refers to single pion production, multiple pion production and inclusive pion production. 



Total cross section 

Differential cross section 

Single 

Isospin inconsistency 
No data at several GeV 
Threshold difficulties 

Isospin needs verification 
Few data mostly below 1 GeV 
Threshold difficulties 
Spin correlation problems 

Multiple 

Similar to Single 

Few data 

Inclusive 

Few data at 1-10 GeV 
Most data at very high energy 

Few data at 1-10 GeV 
Most data at very high energy 
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V. Conclusions 


The present work has considered the availability of pion production data as well as physically motivated, 
parametric fits to such data. It is found that high energy parametrizations do not fit the intermediate energy 
data well, and instead thermal models are more suitable. There are serious gaps in the pion data, with the 
most important being that the data are inconsistent with isospin conservation. More precise data are needed 
as well as inclusive data in the intermediate energy region. Dose contributions from other hadrons also need 
to be evaluated. 
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